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ABSTRACT 

We employ a two-dimensional grid-based hydrodynamic model to simulate 
upper atmospheric dynamics on extrasolar giant planets. The hydrodynamic 
equations of motion are integrated on a rotating, irradiated sphere using a pseu- 
dospectral algorithm. We use a two-frequency two-stream approximation of ra- 
diative transfer to model the temperature forcing. This model is well-suited to 
simulate the dynamics of the atmospheres of planets with high orbital eccen- 
tricity that are subject to widely-varying irradiation conditions. We identify six 
such planets, with eccentricities between e = 0.28 and e = 0.93 and semimajor 
axes ranging from a = 0.0508 A.U. to a = 0.432 A.U., as particularly interesting 
objects for study. For each of these planets, we determine the temperature profile 
and resulting infrared light curves in the 8-/im Spitzer bands. Especially notable 
are the results for HD 80606b, which has the largest eccentricity (e — 0.9321) 
of any known planet, and HAT-P-2b, which transits its parent star, so that its 
physical properties are well-constrained. Despite the variety of orbital param- 
eters, the atmospheric dynamics of these eccentric planets display a number of 
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interesting common properties. In all cases, the atmospheric response is pri- 
marily driven by the intense irradiation at periastron. The resulting expansion 
of heated air produces high-velocity turbulent flow, including long-lived circum- 
polar vortices. Additionally, a superrotating acoustic front develops on some 
planets; the strength of this disturbance depends on both the eccentricity and 
the temperature gradient resulting from uneven heating. The specifics of the re- 
sulting infrared light curves depend strongly on the orbital geometry. We show, 
however, that the variations on HD 80606 b and HAT-P-2b should be readily 
detectable at 4.5 and 8 /im using the Spitzer Space Telescope. Indeed, these 
two objects present the most attractive observational targets of all known high-e 
exoplanets. 

Subject headings: planets: extrasolar — planets: individual (HD 80606b, HD 37605b, 
HD 108147b, HAT-P-2b, HD 118203b, HD 185269b) 



1. Introduction 



In July of 1994, the astronomical community was captivated by the collision of comet 
Shoemaker-Levy 9 and Jupiter. The sudden deposition of energy into the Jovian atmosphere 
caused a dramatic response on a planetary scale, providing remarkable opportunities for 



observation and simulation (Harrington & Deming 2001 Deming & Harrington 2001) and 



fostering the development of a deeper understanding of the structure and behavior of gas- 
giant atmospheres. Within our own solar system, such an occurrence is a one-time event. 
However, the discovery of extrasolar giant planets on orbits with substantial eccentricity 
affords the opportunity to study atmospheric dynamics on worlds where the sudden influx 
of a massive dose of energy occurs on a regular basis. With the advent of Spitzer, direct 
observation of such planets has become feasible. 

The Spitzer telescope has been an unexpectedly useful tool in the study of hot Jupiters, 
and has already been used to make a number of landmark observations of extrasolar planets. 



These began with the detections of the secondary eclipses of HD 209458 b (Deming et al. 



2005) and TrES-1 (Charbonneau et al. 2005), which confirmed that hot Jupiters have low 
albedos and day-side effective temperatures of order T^s ~ HOOK. The TrES-1 observations 
were carried out at 4.5/im and 8/im, and thus provided an estimate of that planet's mid- 
infrared spectral slope. High S/N observations at 16 /im of the secondary transit of HD 



189733 (Deming et al. 2006) produce a similar (1117 ± 42 K) day-side temperature. 



More recently, Harrington et al. (2006) have detected strong orbital phase-dependent 
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brightness variations for v And b at 24 fim, which indicate a strong day-night temperature 
contrast in the far IR, and suggest that heat transfer between the hemispheres is inefficient at 
the level in the atmosphere that generates the observed 24/im emission. Conversely, Cowan et 
al. (2007) found no significant orbital phase variations at 8 fim for 51 Peg b and HD 179949 
b, indicating that for those two planets, surface fiows are effective at redistributing heat 
at the level probed by the 8 fim photosphere. Even more recently, Knutson et al. 2007 
have presented results from a continuous 30-hour Spitzer observational campaign on HD 
189733. The resulting 8/im photometric time-series contains extraordinary (part in 10~^) 
signal-to-noise and clearly shows phase-dependant flux variations that indicate a day-night 
temperature contrast of order 200 K. Of particular interest in the Kuntson et al. (2007) 
time series is a F/ < F >~ 0.0005 increase in planetary flux of ~ 5hr duration just past the 
primary transit. This feature can be interpreted as arising from a hotspot near the dawn 
terminator of the planet. 

In addition to providing infrared photometric time-series data, the Spitzer Space Tele- 
scope can also be used to obtain mid-infrared spectra of extrasolar planets. This has been 
done for the transiting planets HD 209458b (Richardson et al. 2007) and HD 189733b (Grill- 
mair et al. 2007). In both cases, spectra in the 7-14/im range were obtained by subtracting 
the stellar spectrum measured during secondary eclipse from the combined star+planet spec- 
trum obtained just outside eclipse. Model atmospheres of extrasolar planets under strong 
insolation conditions (e.g. Sudarsky, Burrows & Hubeny 2003, Marley et al. 2006) predict 
that H2O absorption features would be prominent within the observed spectral range. Re- 
markably, however, neither planet displayed the predicted signature. This null result can 
be interpreted in several ways. First, it is possible that the atmospheres are dry, perhaps 
because available oxygen is sequestered in CO. This interpretation, however would require 
both planets to have unexpectedly large C to O ratios. Alternately, the effectively featureless 
spectrum could arise from a very clear atmosphere with an isothermal P-T proflle down to 
large optical depth ( Fortney et al. |[2006 ). This interpretation is favored by the fact that the 



hot Jupiters have extremely low albedos (Rowe et al. 2006). As a third possibility, Harring- 



ton et al. (2007) suggest that HD 209458b may be displaying line emission at lO/xm. They 
interpret these features as possibly arising from silicate clouds. Any such clouds, however, 
would need to be highly non-reflective in order to be consistent with the low optical albedo. 

At present, the existing observational results present an incomplete and somewhat con- 
tradictory overall picture. It is not understood how the wind vectors and temperature 
distributions on the observed planets behave as a function of pressure depth, and plane- 
tary longitude and latitude. Most importantly, the effective radiative time constant in the 
atmospheres of short-period planets remains unmeasured, and as a result, dynamical cal- 
culations of the expected planet-wide flow patterns (Cho et al. 2003, Cooper & Showman 
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2005, Burkert et al. 2005, Fortney et al. 2006, Langton & Laughlin 2007, Dobbs-Dixon & 
Lin 2007) have come to no consensus regarding how the surface flow should appear. The 
lack of agreement between the models stems in large part from the paucity of unambiguous 
measurements of the basic thermal structure of the atmospheres of short-period planets. 
Essentially, the problem comes down to the fact that the planets that have been observed 
so far have had their spin periods tidally synchronized to their orbits, so that each spot on 
the substellar surface sees a constant illumination from the parent star. The resulting flow 
pattern will therefore likely be in steady state, and the temperature distribution across the 
planetary surface will depend on both the radiative timescale as well as the windspeeds at 
optical depth unity. 

In this paper, we discuss our modeling program for determining the hydrodynamic 
flows on the surfaces of planets traveling on significantly eccentric orbits. Such planets 
are not spin-synchronous. Tidal effects cause the rotation rate to approximate the orbital 
angular velocity at periastron; the rotation is therefore pseudo- synchronous and the spin 



period is a sensitive function of the eccentricity (Hut 1981). In contrast to the typical spin- 
synchronous hot Jupiter, on which one hemisphere is perpetually illuminated while the other 
receives no insolation, the entire surface of an eccentric planet receives radiation. Pseudo- 
synchronization does, however ensure that a single hemisphere receives the preponderance of 
the illumination during the periastron passage. Additionally, the intensity of the radiation 
received varies strongly over the course of an orbit. The insolation of planets on eccentric 
orbits therefore exhibits significant time variation. 

To our knowledge, other modeling efforts to date have focused exclusively on exoplanets 
that are in tidally synchronized orbits with a permanent day-night side asymmetry and 
essentially steady state irradiation conditions. With such bodies, it is difficult to disentangle 
the contributing effects of the effective radiative time constant in the planetary atmosphere, 
and the hydrodynamic flow velocities on the surface. In this paper, we propose that these 
ambiguities can be resolved by observing planets that have eccentric orbits, and which hence 
show different hemispheres (or partial hemispheres) to the star, and which experience a 
highly variable total insolation. We identify the highly eccentric planet HD 80606 b and 



the transiting planet HAT-P-2 b (Bakos et al. 2007) as being a prime targets for such 



observations. We propose that by monitoring these planets during their periastron passages, 
we can improve our overall understanding of the atmospheric properties of hot Jupiters. 
The planet-wide response to the rapid increase in stellar radiation will be readily observable 
in the 4.5/im and 8/im bands, and will provide a direct simultaneous measurement of the 
radiative time constant at two separate levels in the atmosphere. 

The plan for this paper is as follows. In §2, we describe our numerical model as applied 
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to extrasolar planets. In §3, we discuss the characteristics of the six planets considered in this 
paper, as well as relevant stellar properties of their parent stars. In §4, we present the results 
of this model as applied to these planets. In §5, we discuss the observational consequences of 
these results, in particular concentrating on the infrared light curves in the 8 micron Spitzer 
band. We show that in general, the periastron passage produces a distinct flux peak, and 
that the shape of the overall curve can be used to constrain planetary characteristics such 
as rotation rate and atmospheric radiative time constant. We conclude in §6. 



2. Numerical Model 

We begin with the primitive hydrodynamic equations of motion in two dimensions, 
assuming hydrostatic equilibrium: 



-v-Vp-pV-v (1) 



dp 
di 

dv 1 

m = (2) 

We use the equation of state for an ideal gas: p = (7 — l)e = pRT, where R is the specific gas 
constant and the adiabatic constant 7 = 1.4 for a diatomic gas. (We actually use a value of 
7 = 1.389, appropriate for molecular hydrogen at 900 K.) Setting the specific heat capacity 
to Cv = R/{l — 1), we can substitute 

T=~ (4) 

Oy p 

into equation ([3]) to obtain the hydrodynamic equations in terms of the dynamical variables 
T, V, and p: 



V ■ VT - (7 - 1)TV ■ V (5) 



dT 
~dt 

^ = _v . Vv - ^Vp - i?VT (6) 



dp 
di 



-V ■ Vp - pV ■ V. (7) 



In order to obtain a model appropriate to planetary atmospheric flows, we add a thermal 
forcing term to ^ and a Coriolis term to (|6]), obtaining the equations of hydrodynamic 
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motion on an irradiated, rotating sphere: 

(7 - 1)TV ■ V + frad (8) 

— Vp - RVT - 2nrot sin e{n x v) (9) 

pV ■ V, (10) 

where frad gives the thermal forcing, Qrot is the angular frequency of the planet's rotation, 
6 is the latitude, and n is a unit vector normal to the planet's surface. 

We calculate frad using a one-layer, two-frequency radiative transfer scheme. For hot 
Jupiters, the bulk of the energy input occurs near the peak of the stellar spectrum in the 
visible Amax ~ 2.73/T^ ~ O.Sp. Energy balance, on the other hand, is maintained via thermal 
re-radiation in the infrared Amax ~ 2.73/Tpi ~ 3fi. The separation of the corresponding 
blackbody distributions motivates our use of a two-frequency scheme. We take the mean 
opacity for incident stellar radiation to be ki, while the mean opacity for outgoing long- wave 
radiation is /c2- In this greatly simplified model, a thin layer of the atmosphere with thickness 
dz and at a pressure depth p will absorb a flux equal to 



dT 

~dt 
(9v 

dp 
di 



-V- VT 
-V ■ Vv 

-V ■ Vp - 



dF=( hiM + k2aX^ - harA p dz, (11) 
\ cos a / 

where F^,{p) is the incident stellar flux at a pressure depth p, a is the zenith angle of the 

star, Tn is the temperature of the layer being modeled in the absence of solar heating, 

and a = 5.67 ■ 10~^ W m~^ is the Stefan-Boltzman constant. T„ is the result of the 

combination of the intrinsic planetary flux due to tidal heating and the heating supplied by 

the atmosphere at high optical depth, and is calculated using a one-dimensional two-stream 



approximation to radiative transfer similar to that used in Iro et al. (2005). The details of 
the Tn calculation are discussed in greater detail in §3. 

Energy conservation requires that dF = pCpfraddz, so the thermal forcing must be 

/...= (^)(^ + Mr„'-fe.T'). (12) 

\CpJ \ COS a / 

We may approximate -F*(p) by ignoring the variation in ki with temperature and pres- 
sure; this gives 

F,(p) = F,(0)e"^. (13) 



But F^{0) is simply the flux at the top of the atmosphere, which is 



F,(0) = (1 - A) ( ) cosa, (14) 
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where A is the Bond Albedo, L^, is the stellar luminosity, and a = a{t) is the time-dependent 
distance between the star and the planet. Then the thermal forcing is 

frad = (a:i(1 - A) + k,aT,', - k,aT'^ . (15) 

With a little bit of algebraic manipulation, this can be placed in a more transparent form: 

frad = HTe,-T'), (16) 

where /5 = <7k2/Cp and the equilibrium temperature T^g follows 

with (tT^^ = (1 — A)L^,/{4:TTa'^) and x = exp{—kip/ g). The forcing that results from this 
scheme strongly depends on the choice of ki, k2, and p. Motivated by many-layer radiative 



models of HD 209458 b ( |Iro et al.||2005| , we choose ki = 2 ■ 10'^ m^ kg-\ A;2 = 4 ■ 10"^ 
m^ kg~^, and p = 250((7/10 m s~^) mbar, where g is the acceleration due to gravity at the 
planet's surface. 

While we believe this model represents an improvement over earlier work, it is never- 
theless unable to account for a number of potentially important effects. Two-dimensional 
hydrodynamics obviously cannot account for vertical atmospheric flows. Since the infrared 
photosphere is well inside the radiative zone for the planets under consideration, it is likely 
that convective effects are negligible. However, the inability of the model to account for 
vertical expansion under heating could cause some error. 

In our radiative model, we make a number of simplifications are likely to cause dis- 
crepancy between our predictions and observations. We ignore the variation of opacities 
with temperature and pressure, which is known to be significant. We furthermore do not 
consider the variation in opacity with wavelength within the infrared regime; this implies 
that, for example, the 24- micron photosphere and the 3- micron photosphere occur at the 
same pressure level. The planets in our simulations experience large temperature variations. 



which may lead to the formation and dispersal of different cloud species (Sudarsky, Burrows 



& Hubeny 2003 Marley et al. 2007). We do not, however, include the formation of clouds 



in our model. 

Despite these shortcomings, this model represents a significant improvement over earlier 
simulations in a number of areas. Shallow water models can achieve high resolution, but 
also have difficulty modeling certain qualitative aspects of the flow; temperature waves redis- 
tribute heat far more efficiently than is realistic. Since our code employs fully compressible 
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hydrodynamics, we are able to simulate heat transfer within atmospheric flows with greater 
realism than is possible using incompressible shallow-water dynamics. 

Furthermore, our radiative scheme provides a more realistic model of the thermal forc- 



ing than the Newtonian relaxation employed in earlier hydrodynamic simulations (Cho et al. 



2003 Cooper & Showman 2005 Langton & Laughlin 2007). While Newtonian relaxation is 



a reasonable approximation when the temperature is close to the equilibrium temperature, 
it becomes dramatically less accurate as the temperature perturbations grow larger. In the 
forcing regimes expected on extrasolar planets, the Newtonian approximation can overesti- 
mate the actual rate of temperature change by as much as a factor of 3. Our radiative model 
therefore allows us to treat the stellar heating more accurately than previous simulations; 
this is particularly important in cases where the thermal forcing varies significantly with 
time, either due to non-synchronous rotation or to a highly eccentric orbit. 



Equations (\8v through (10) are numerically integrated using the vector spherical har- 



monic transform method (Adams & Swarztrauber 1999). We use a longitude-latitude grid 



with T171 resolution (512 x 257), sufficiently fine to resolve the details of turbulent fiow. Nu- 
merical stability is maintained using fourth-order hyperviscosity ( Polvani et al.|[2004 ), which 
prevents the spurious addition of energy at small scales. This is accomplished by adding 
a term of the form B'V^(j) to each of the time-evolution equations, where (p is the evolving 
quantity. To minimize the effect of this numerical dissipation on the flow evolution, B is 
generally chosen to be as small as possible, while still preventing unphysical spectral ringing 



or numerical instability. In general, this requires B 



where a is the planetary 



radius. This results in a decrease in of ~ 5 % per time step at the smallest spectral scales; 
with such strongly-varying heating, anything less allows pernicious behavior. 

The planetary atmosphere is initially at constant pressure, isothermal, and nearly at 
rest, with only a small random initial flow to break the unphysical north-south symmetry. 
Since the resulting flow is often turbulent, the speciflc evolution of the dynamical variables 
changes from one run to the next. Nevertheless, the qualitative properties of the flow are not 
altered by small changes to the initial wind flow. The initial temperature is chosen so that it 
approximately equals the average equilibrium temperature; this prevents the initial thermal 
forcing from inducing the large winds that result if, for example, the initial temperature is T„ 



(Cooper & Showman 2005). While other groups have obtained supersonic flow, these seem 



to occur in simulations where either the night-side equilibrium temperature is unusually low 
( Dobbs-Dixon fc Lin|2007 ) or the planet is initially cold ( Cooper &: Showman|2005 ) . The flrst 
situation maintains an artiflcially high equilibrium temperature gradient; the second greatly 
amplifles the effect of the initial heating, particularly if the stellar heating is not introduced 
gradually. In either case, there is good reason to believe that the actual wind speeds may be 
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subsonic. The planet is allowed to evolve from its initial state until a dynamical equilibrium 
is reached - typically ~ 100 days - and then the temperature, density, wind velocity, and 
relative vorticity are recorded. 



3. Planetary Characteristics 



Our code requires a number of planetary characteristics as input parameters. The 
mass and radius of the planet must be specified in order to determine the surface gravity 
and therefore the appropriate pressure at the infrared photosphere. The orbital period, 
semi-major axis, and eccentricity of the planet's orbit are required to calculate the time- 
dependent irradiation conditions. Additionally, the longitude of periastron is needed to 
produce the light curves. With the exception of the planetary radius, all of this information 



can be found in Butler et al. (2006). The radii are determined using the planetary structure 



model of Bodenheimer et al. (2003). To determine correctly the mass and the light curve, it 



is also necessary to know the orbital inclination. In general, however, this can be constrained 
only in the case of a transit. Nevertheless, inclinations close to 90° are more probable. We 
therefore take the minimum value of the planetary masses; that is, sini = 1. The effect 
of the planetary mass enters the simulation primarily by determining the magnitude of the 
gravitational acceleration at the planetary surface. Test runs of our model varying gravity 
by up to a factor of 10 produce little change in the resulting flow patterns. Since a two- 
dimensional model of a radiative atmosphere necessarily assumes hydrostatic balance and 
vertical stratification, the small impact of gravity is not surprising. In any case, it is clear 
that in our model, the actual mass of the planet is of negligible importance; the minimum 
mass assumption therefore introduces no significant inaccuracies. 

The planet's rotation rate is also an important parameter, necessary to determine both 
the substellar point and the strength of the Coriolis effect. In the case of a close-in cir- 
cular orbit, planetary rotation is expected to be synchronous, so that the orbital period 
and the rotation period are identical. If the orbit is even slightly eccentric, however, spin- 
synchronization is broken. Tidal effects are strongest at periastron, and therefore the pre- 
ferred rotation rate is approximately the angular frequency at periastron. The result is a 
state of spin pseudo-synchronization, wherein the rotation rate is a well-defined function of 



the eccentricity (Hut 1981): 



' rot 



^orb 




Recently, Ivanov & Papaloizou (2007) have advanced an alternate theory for spin- 



- 10 - 



synchronization which predicts that planets in highly eccentric orbits will have spin angular 
velocities that are close to 1.55 times the circular orbit angular velocity at periastron. As 
will be clear from our simulations in this paper, long-duration observations with the Spitzer 
Space telescope can in principle be used to distinguish between the two models for spin 
pseudo-synchronization. In this initial work, however, we adopt the Hut (1981) theory for 
our assumed planetary spin rates. 

In addition to the planet's physical and orbital characteristics, we also must determine 
the equilibrium night-side temperature, T„ at the simulation pressure depth p = 250(^^/10 
m s~^) mbar. This represents a combination of the planet's intrinsic luminosity due to tidal 
dissipation and gravitational contraction and the redistribution of incident stellar radiation 
at high optical depths. With the exception of HD 80606 b, which is characterized by an 
extremely high tidal luminosity and relatively low average insolation, the night side temper- 
atures of the planets here considered are largely determined by the incident radiation. 

We determine T„ using a one-dimensional plane-parallel two-stream algorithm with = 
100 vertical zones. In keeping with the radiative model described in §2, we consider only two 
frequencies and neglect the variation in opacities with temperature and pressure. The results 
of these calculations are therefore not precise, but are likely to give reasonable estimates of 
the atmospheric temperature in the absence of direct stellar heating. At the upper boundary, 
we impose a downward flux equal to the time-averaged stellar radiation 

^^^^^1Q^{^)^1M^/T^' ^^^^ 
At the lower boundary, which we set at p = 5.0 ((7/ 10m s~^) bar, we require a net upward flux 
consistent with the planet's intrinsic luminosity. We then allow the atmosphere to come to 
radiative equilibrium, thereby determining the temperature at the base of the atmosphere. 
We then maintain the lower boundary at this temperature, and recalculate the equilibrium 
temperature at each layer in the absence of stellar radiation. This quantifies the response 
of the upper atmosphere to heating by high-temperature interior layers. We take T„ to be 
the atmospheric temperature at the layer corresponding to the pressure depth used in our 
simulations. The planetary parameters, including T„ and Trot, are summarized in Table [1} 

Finally, the characteristics of the parent star are necessary to determine the magnitude 
of the incident stellar radiation, as well as the relative planet-star infrared flux as a function 



of wavelength. From Butler et al. (2006), we obtain the stellar mass, radius, temperature, 
and luminosity, although our model only requires the temperature and luminosity. These 
data are compiled in Table [2] 

We wish to focus particular attention on two planets: HD 80606 b and HAT-P-2 b. 
With e > 0.9, HD 80606 b is the most eccentric planet known by a significant margin. 
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Because both the rotation rate and the characteristics of the periastron passage depend 
strongly on the eccentricity, it is especially important to obtain an accurate value. Using 



both the 55 radial velocities published in support of the initial detection (Naef et al. 2001) 



and a further 46 radial velocities published by Butler et al. (2006), we obtain best-fit values 



for the relevant planetary characteristics: the orbital period Tort = 111.4298 ± 0.0032d, the 
eccentricity e = 0.9321 ± 0.0023, and the longitude of periastron w = 300.264° ± 0.35°. 

HD 80606's mass is estimated to lie between 0.9 and 1.1 Mq (Naef et al. 2001, Butler 
et al. 2006). Adopting M = IMq, we derive Mpisin(z) = 4.18 Mjup. Our planetary models 
(Laughlin et al. 2005) predict a radius of 1.10 ± 0.02i?jup. The uncertainty in Rp\ is small 
because the planet lies in the mass regime where the radius is an extremely weak function 
of the mass. Furthermore, the planet is sufficiently massive enough that while it experiences 
considerable internal heating associated with its strong tidal dissipation, the heating (which 
we estimate below to be L ~ 10^^ ergs) is expected to have a negligible impact on the planet's 
physical radius (Gu, Bodenheimer, & Lin 2004). The well-determined radius is important 
because it largely eliminates an unknown free parameter when we model the hydrodynamic 
response. 

Because of the extreme eccentricity and relatively large semi-major axis of HD 80606 b's 
orbit, tidal heating becomes a significant contribution to the overall planetary temperature. 



Wu & Murray (2003) posit that HD 80606 b's large eccentricity was initially produced by 
resonant interaction with HD 80607, a binary companion to HD 80606 with a wide separation 
~ 1000 A.U. Wu and Murray show that if HD 80606 b initially formed in an ainit = 5 A.U. 



orbit that was highly inclined to the binary plane, then the Kozai resonance (Kozai 1962) 
would have quickly driven the planet's eccentricity to e ~ 0.993. In the Wu-Murray theory, 
successive cycles of high eccentricity and attendant tidal orbital decay acted to decrease the 
semi-major axis. Once a dropped below ~ 3 A.U., the Kozai resonance was destroyed by 
general relativistic precession. Over the last few Gyr, the planet has undergone continuous 
tidal dissipation, which has decreased a and e to their current values. By considering the 
amount of orbital energy lost over this time period, [Wu &: Murray] arrive at a tidal luminosity 



Ltid 10^^ erg. For the sake of comparison, the average incident stellar radiation absorbed 
by the planet is Lrad = 3 ■ 10^'' ergs. It is likely, therefore, that the dominant contribution 
to Tn on HD 80606 b arises from eccentricity tides rather than convective redistribution of 
incident radiation. 

HAT-P-2 b also deserves special consideration. It was discovered via the transit method 
by Bakos et al. (2007). Because it transits its parent star, its mass M = 8.17 Mj it 0.72Mj 



and radius R = l.lSRj it O.lQRj are well-constrained. In particular, the sini degeneracy is 
broken, since the inclination i = 90.0° ± 1° has been directly measured. Furthermore, we 
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have direct determination of the radius; this is particularly valuable given the divergence 
of observed radii from models in a number of cases ( [Sato et al. 2005 Bodenheimer et al. 



2003t [Laughlin et al.|[2005| . HAT-P-2b orbits a bright (V=8.7) star with a period of 5.63 d; 
furthermore, its orbital geometry [w = 185°), ensures that the periastron passage is midway 
between the primary and secondary transits and that all three events are observable within 
a 36-hour period. HAT-P-2b is therefore uniquely well-suited for observation. 



4. Dynamical Evolution of the Planetary Atmospheres 

Representative atmospheric temperature distributions for each planet considered in this 
paper, with the exception of HD 37605 b, are presented in Figure [T} Our simulations demon- 
strate that eccentric exoplanets exhibit great diversity in their atmospheric flow patterns and 
temperature distributions, with eccentricity, periastron distance, orbital period, and stellar 
luminosity all playing important roles. Nevertheless, all of the planets here considered dis- 
play significant subsolar heating at periastron, with temperature increases ranging from 300 
K to 700 K. Planets with high eccentricities and strong periastron heating may also develop 
high-amplitude superrotating acoustic waves in response to a high temperature gradient near 
the eastern terminator. Due to the temperature dependence of the sound speed, the acoustic 
waves emitted inside hotter regions overtake those emitted at the cooler terminator, causing 
a region of very high density only a few zones wide. This "acoustic compression" produces a 
front propogating eastward, maintaining its intensity as it is repeatedly focused by reflections 
off of the cooler, denser polar air. 

It is not immediately obvious why our model predicts the formation of these acoustic 
disturbances on some planets but not on others. The orbit of HD 37605 b is considerably 
more eccentric (e = 0.737) than that of HAT-P-2 b (e = 0.507), yet HAT-P-2 b exhibits 
a clearly defined acoustic disturbance which is absent from HD 37605 b. The periastron 
flux received by HD 80606 b is smaller than that received by HD 118203 b and HD 185269 
b, which orbit more luminous stars. Furthermore, HD 80606 b has a comparatively small 
temperature variation and correspondingly slower wind speeds. Yet the stellar heating re- 
ceived by HD 80606 b generates large-amplitude acoustic waves, whereas the flow on HD 
118203 b and HD 185269 is primarily vortical, with no large-scale acoustic fronts. While we 
have not investigated this issue in detail, our calculation of predicted light curves suggests 
that the subsolar heating rate at periastron is the determinating factor in the formation of 
a strong acoustic front. Sufficiently rapid heating at periastron produces a large local tem- 
perature gradient which, as described above, generates a region of high density and pressure 
which then propagates acoustically; slower periastron heating causes a shallower tempera- 
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ture gradient, less acoustic compression, and therefore any acoustic waves are much lower 
amplitude. 

In order to elucidate the dynamics of the flow itself, we show in Figure [2] the vorticity 
distribution on five of the six planets modeled; again, we omit the results for HD 37605 b. 
The distributions are shown at the same time and, in the case of the globe view, from the 
same vantage point as in Figure [TJ A quiver plot superimposed on the equirectangular map 
shows the direction and relative wind speed of the flow itself. 

From both the temperature and vorticity distributions, it is clear that with the excep- 
tion of HD 80606 b, each planet exhibits large-scale turbulent eddies. Given the extreme 
eccentricity of this planet's orbit, this is perhaps surprising. However, the lack of turbulence 
can be explained by the relatively low temperature contrast (which we take to be Tmax — Tmin, 
where Tmax and Tmin are respectively the global maximum and minimum temperatures at a 
given moment) seen on HD 80606 b. While this temperature contrast reaches 500 K momen- 
tarily, the rapid periastron passage prevents the planet from receiving enough heat to sustain 
these temperature differences; within 72 hours of periastron, the temperature difference has 
fallen to 250 K, continuing to drop after that. In comparison, the temperature difference on 
HAT-P-2b exceeds 1100 K at periastron, with an average value of 540 K. While we have not 
run a sufficient number of simulations to be precise, it appears that sustained temperature 
differences greater than 250 to 300 K are required to generate large-scale vortical struc- 
ture. It is important that this temperature contrast be prolonged; while brief exposure to a 
large temperature gradient induces divergent flow, it appears that it takes some time for the 
Coriolis force to produce a large and persistent rotational component. 

It is interesting to note here the significance of T„. For a particular amount of stellar 
heating, as T„ increases, the induced temperature gradient becomes smaller, and both wind 
speed and turbulence decrease. We find, for example, that a lower value of tidal heating on 
HD 80606 b, leading to T„ = 400 K provides a sufficiently cool background that HD 80606 
b develops a temperature difference which remains over 300 K for 10 d after periastron, 
thereby producing turbulent flow comparable to the other planets. Similiarly, a test run on 
HD 118203 b with T„ = 1960 K (a change with absolutely no physical justification!) reduced 
the maximum temperature contrast on that planet to 200 K, stifling the development of the 
hemisphere-scale vortices so prominent in our earlier simulation. 

The intrinsic temperature of the planet then rises to a question of first-order importance. 
While the infrared observations resulting from the evolution of a particular temperature dis- 
tribution are not sensitive to small-scale differences in the flow, a global-scale qualitative 
shift from laminar to turbulent flow fundamentally alters the temperature distribution, with 
possibly significant observational impact. Given the importance of tidal heating in determin- 
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ing the intrinsic luminosity, in addition to the unsuitabihty of an equihbrium model of tides 
to predicting the tidal dissipation on these highly dynamic objects, further development of 
a dynamical theory of tides seems to be a promising field for future investigation. 



5. Infrared Emission Curves 

The flux F\ emitted by a planet with some temperature distribution T(0, 9) between 
wavelengths A and A + c/A is 

Fx = j BxiT) cos a cos OrldedcpdX, (20) 

where Bx{T) is the Planck function for emission at a wavelength of A by a blackbody at 
temperature T, and a is the angle between the surface normal and a vector directed from 
the planet to the observer. The integral is performed over the hemisphere visible to the 
observer, that is, where cos a > 0. In the case of a grid-based model, this can be discretized 
as follows: 

Fx = Y^ Bx{T) cos a cos OrpAOA^, (21) 

where 

Bx{T) = [ Bx{T)dX. (22) 



We apply equation (21) to the temperature distributions considered in §4 to produce plane- 
tary infrared light curves; the results are shown in Table [3] and Figures |3||4| and|5j In Figures 
|4] and |5j the small points on the orbital diagram show the position of the planet at one-hour 
intervals; the larger globes (labeled 1 through 4) show the position of the planet at times 
t = — 19h, Oh, +19h, +38h relative to periastron. The black poles on these globes show the 
rotation of a stationary point on the planetary surface; the gray poles indicate the amount of 
rotation that has occurred over a 19-hour period. The hemispheric plots in the right portion 
of the figures show the atmospheric temperature distribution at the corresponding times. 
The light curves shows the planets' infrared emission in the 8 /xm band from 25 h before to 
50 h after periastron. 

These light curves vary widely in their shape, range, and time scales, as is expected 
given the range of irradiation conditions and orbital parameters. All show some degree of 
interesting structure; none can be well-fit by the simple sinusoidal shape predicted by earlier 



models (Cooper & Showman 2005; Langton & Laughlin 2007). We note two particularly 



promising candidates for observation. HD 80606 b, which reaches a peak flux F/F* = 8.0 ■ 
10~^ in the S-fim band, will be observed by the Spitzer Space Telescope during its periastron 
passage in November 2007. HAT-P-2 b, which exhibits even larger flux variation over the 
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course of its orbit, provides an equally attractive target for observation. The flux variations 



for both objects are sufficiently large to be comfortably resolved by Spitzer (Charbonneau 



et al. 2005 Harrington et al. 2006 2007). The other four planets, HD 37605 b, HD 108147 



b, HD 118203 b, and HD 185269 b, present significantly greater observational challenges. 
However, it is not infeasible that the light curve variations on some of the brighter planets - 
HD 118203 b and HD 108147 b - could be resolved with sufficiently long integration times. 



6. Conclusion 

The advent of the Spitzer Space Telescope and the discovery of short-period planets 
on significantly eccentric orbits has presented astronomers with the opportunity to observe 
time-dependent weather patterns on extrasolar planets. Giant planets on eccentric orbits 
that experience strong variations in insolation have the potential to tell us a great deal about 
the structure, dynamics and chemistry of Jovian-like atmospheres. Eccentric planets can be 
profitably compared with short period planets on spin-synchronous circular orbits that are 
likely to have largely steady-state patterns of surface fiow, and which are currently being 
observed with Spitzer. 

Furthermore, when eccentricities are extreme - as in the case of HD 80606 b - we are 
afforded the ability to directly observe planets that spend the bulk of their orbit in a region 
where the radiation environment is similar to that experienced on Earth. With current 
instruments, such planets would be completely beyond the possibility of direct observation. 

Eccentric planets also have the best potential to give us an understanding of the dissi- 
pative mechanisms at work in planetary interiors. A planet such as HAT-P-2 b should be 
generating a very significant luminosity as a result of the strongly time-dependent tidal force 
on the planet. By observing the depth of the primary transit in the infrared, one can get a 
direct measurement of the planetary effective temperature, which can then be used to infer 
the amount of dissipation within the planet, and to put interesting constraints on the tidal 
quality factor, Q associated with dynamical (as opposed to equilibrium) tides. 

To take a first step toward understanding what to expect from the infrared light curves 
of short-period eccentric planets, we have created a two-dimensional hydrodynamic model 
with radiative forcing appropriate to the time-varying illumination received for a specified 
orbit. Our model handles radiation transfer in an approximate, but hopefully qualitatively 
correct manner by separating the transport of optical and IR radiation. We apply this 
model to six known giant planets on eccentric orbits with e > 0.3. The irradiation intensity 
on these planets varies by a factor of three for the least eccentric planets up to a factor 
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of 800 for HD 80606b. The effect of this tremendous variabihty on atmospheric flows is 
profound; the dominant feature in the evolution of each atmosphere stiidicd is the subsolar 
heating through the perastron passage, with the subsequent flow being primarily driven by 
this dramatic forcing. For each planet, we find significant spatial and temporal variability 
in the state of the atmosphere, with turbulent fiow being a common feature. Additionally, 
rapid increase in the thermal forcing leads to superrotating acoustic shocks, as seen in the 
atmospheres of HD 80606 b and HAT-P-2b. 

Among the planets that we have modeled, HAT-P-2b and HD 80606b stand as the most 
promising candidates for time-series infrared observations in Spitzer's S-fira band, and among 
these two planets, HAT-P-2 b is the best candidate. HAT-P-2b's secondary transit provides 
an opportunity for an overall flux calibration, and its orbital geometry (with zu — 185°) 
places the most interesting features of the light curve on a relatively short 35-hour segment 
within the full 5.6-day orbit. During this time span, we predict that the relative infrared 
flux will vary from 3.2 x lO"'' to 9.0 x lO"''. 

HD 80606b, on the other hand, is predicted to produce a relative infrared flux ranging 
from 3.8 x 10~^ to 8.0 x 10~^. This variation easily satisfies Spitzer's sensitivity hmit, and 
for this planet, the observational situation is greatly aided by HD 80607, the wide binary 

companion of near-identical spectral type to the parent star HD 80606. Billions of years ago, 
HD 80607 enabled Kozai migration. Today, it provides a stable photometric reference to 
gauge the rapid infrared fiux rise from the planet during periastron passage. A measurement 
of the slope of this rise will allow a direct estimate of the atmospheric radiative time constant 
at the A.5-/J and 8-/x photospheres. 

A short-period Jovian planet on an eccentric orbit likely presents one of the Galaxy's 
most thrilling sights. One can imagine, for example, how HD 80606 b appears during the 
interval surrounding its hair-rising encounter with its parent star. The blast of periastron 
heating drives global shock waves that reverberate several times around the globe. From 
Earth's line of sight, the hours and days following periastron are characterized by a gradually 
dimming crescent of reflected starlight, accompanied by planet- wide vortical storms that fade 
like swirling embers as the planet recedes from the star. It's remarkable that we now have the 
ability to watch this scene (albeit at one-pixel and two-frequency resolution) from a vantage 
several hundred light years away. 

Indeed, the observational and theoretical allure of short-period non-synchronous planets 
is worthy of much more detailed study than we have presented here. Our next steps on this 
problem, therefore, will be to (1) upgrade our radiative transfer scheme to allow many 
more frequency intervals to be modeled, and (2) incorporate the higher-resolution radiative 
transfer into a full three-dimensional hydrodynamic treatment. 
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Torb{(i) 


TroL{d) 


a (AU) 


e 


i 


zu 


Alpi sin t 
Mj 


Rp, 
Rj 


Tn (K) 


HD 80606 b 


111.449 


1.72 


0.432 


0.9321 




300.3° 


4.18 


1.1 


720 


HD 37605 b 


54.22 


6.76 


0.261 


0.737 




211.6° 


2.86 


1.1 


320 


HD 108147 b 


10.90 


3.50 


0.102 


0.53 




76.7 


0.261 


0.9 


670 


HAT-P-2 b 


5.63 


1.96 


0.0685 


0.507 


90° 


185° 


8.17 


1.18 


1010 


HD 118203 b 


6.134 


3.85 


0.0703 


0.309 




155.7 


2.14 


1.1 


960 


HD 185269 b 


6.8380 


4.61 


0.0766 


0.280 




100° 


0.91 


1.0 


935 



Table 1: Characteristics of planets studied. Torb and Trot arc the orbital and rotational periods, 
respectively. T„ is the equilibrium night-side temperature (see §2 and 3). The other variables 
have their usual meanings. The inclination i is known only for the transiting planet HAT-P-2 
b. 







R*/ Rq 






V 


HD 80606 


1.0 


1.05 


5573 


0.947 


9.06 


HD 37605 


0.800 


0.967 


5391 


0.707 


8.67 


HD 108147 


1.19 


1.29 


6156 


2.14 


6.99 


HAT-P-2 


1.35 


1.80 


6290 


4.25 


8.71 


HD 118203 


1.23 


2.13 


5600 


4.00 


8.05 


HD 185269 


1.28 


2.01 


5980 


4.61 


6.67 



Table 2: Characteristics of parent stars. Only the radius i?*, the temperature T*, and the 
luminosity arc directly used in computing the light curves; however, the mass M* and 
the apparent magnitude V are also included for reference. 





Tfjiin (K) 






■^max / 


AP/P 


v^ax (m s--^) 


Vrnis (m S--^) 


HD 80606 b 


700 


1200 


3.8 • 10-4 


8.0 • 10-4 


0.2 


670 


200 


HD 37605 b 


340 


690 


7.1 • 10-5 


2.6- 10-4 




300 


70 


HD 108147 b 


680 


1400 


1.7-10-4 


4.0 • 10-4 




1060 


370 


HAT-P-2 b 


950 


2170 


3.2 • 10-4 


9.0 • 10-4 


0.5 


1460 


360 


HD 118203 b 


920 


1770 


2.1 • 10-4 


4.5 • 10-4 




1140 


390 


HD 185269 b 


900 


1760 


1.8- 10-4 


2.9- 10-4 




720 


1710 



Table 3: Summary of simulation results. T„jj„ and Tmax arc the global minimum and max- 
imum temperatures, respectively; Fmin/F* and F„iax/F* are the minimum and maximum 
relative planetary fluxes in the 8 micron band. AP/P gives the intensity of the over- pressure 
front on planets where an acoustic disturbance is present. The maximum and root-mean- 
square wind speeds are shown in the last two columns. The rms wind speed for HD 80606 b 
varies significantly over the course of the simulation due to the high eccentricity of its orbit; 
the peak value of Vrms is given for this planet. 
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Fig. 1. — Temperature distributions at representative (yet photogenic) positions in the orbit 
of each planet modeled. The globes on the left show a three-dimensional representation 
of the temperature; the plots on the right show an equirectangular representation. The 
temperature ranges in each plot follows: HAT-P-2 b: 950 K to 1500 K; HD 80606 

b: 710 K to 890 K; HD 185269 b: 930 K to 1530 K; HD 108147 b: 700 K to 1390 K; HD 
118203 b: 970 K to 1770 K. 
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Fig. 2. — Relative vorticity distributions and wind flow patterns for each planet modeled. 
The globes on the left show a three-dimensional representation of the (absolute value of) 
relative vorticity. The plots on the right show an equirectangular plot of the wind velocity 
field superimposed on the vorticity distribution. The flow pattern is taken at the same time 
as the temperature distribution in Figure [T| The maximum wind speeds in each plot are 
as follows: HAT-P-2 b: 1050 m s"^ HD 80606 b: 380 m s"^; HD 185269 b: 1410 m s'^ 
HD 108147 b: 930 m s~^; HD 118203 b: 1100 m s~^. The maximum relative vorticities are: 
HAT-P-2 b: 1.1-10-^ s'^ HD 80606 b: 1.8-10-^ s"^; HD 185269 b: 1.2-10-^ s"^ HD 108147 
b: 1.3-10~^ s"^] HD 118203 b: Ll-lO""^ s"^ (Note that the relative vorticity at a particular 
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Fig. 3. — 8-/im infrared emissions for all six planets modeled over a 250-hour period. The 
cyan curve shows the light curve for HAT-P-2 b; the black curve, HD 80606 b; the red curve, 
HD 118203 b; the purple curve, HD 108147; the cerulean curve, HD 185269 b; the navy blue 
curve, HD 37605. These light curves were prepared under the assumption that i = 90°. 
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Fig. 4. — A detailed representation of the results of the HAT-P-2 b simulation. The 
small points on the orbital diagram show the position of the planet at one-hour inter- 
vals; the larger globes (labeled 1 through 4) show the position of the planet at times 
t = — 19h, Oh, +19h, +38h relative to periastron. The black poles on these globes show the 
rotation of a stationary point on the planetary surface; the gray poles indicate the amount of 
rotation that has occurred over a 19-hour period. The hemispheric plots on the right show 
the atmospheric temperature distribution at the corresponding times. The light curve shows 
the planet's infrared emission in the 8 fim band from 25 h before to 50 h after periastron. 
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Fig. 5. — A detailed representation of the results of the HD 80606 b simulation. The 
small points on the orbital diagram show the position of the planet at one-hour inter- 
vals; the larger globes (labeled 1 through 4) show the position of the planet at times 
t = — 19h, Oh, +19h, +38h relative to periastron. The black poles on these globes show the 
rotation of a stationary point on the planetary surface; the gray poles indicate the amount of 
rotation that has occurred over a 19-hour period. The hemispheric plots on the right show 
the atmospheric temperature distribution at the corresponding times. The light curves show 
the planet's infrared emission in the 8 /im band from 25 h before to 50 h after periastron. 
The red curve corresponds to an orbital inclination of i = 90° (in which the orbit is being 
viewed edge-on), the purple to z = 60°, the blue to i = 30°, and the black to z = 0°. 



